Tensor-entanglement renormalization group approach to topological phases 
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The tensor-entanglement renormalization group approach is applied to Hamiltonians that realize 
a class of topologically ordered states - string-net condensed states. We analyze phase transitions 
between phases with and without string-net condensation. These phase transitions change topo- 
logical order without changing any symmetries. This demonstrates that the tensor-entanglement 
renormalization group approach can be used to study the phase diagram of a quantum system with 
topologically ordered phases. 



Introduction: A mean field approach could potentially 
be very useful for understanding and analyzing topo- 
logical phases[l] and topological phase transitions. [2-8] 
The main challenges in developing such an approach are 
(a) finding a class of mean field states that can describe 
topological phases and (b) finding a simple way to cal- 
culate the physical properties, such as average energy, 
of these states. The first problem can be solved with 
a general class of trial wave functions known as "ten- 
sor product states" (or alternatively "projected entan- 
gled pair states"). [9, 10] Indeed, one can show that ten- 
sor product states (TPS) can describe all the string-net 
condensed states constructed in [11], and hence all non- 
chiral topological phases. [12] The second problem can 
be solved with the tensor-entanglement renormalization 
group (TERG) method[13] which allows one to calculate 
correlation functions (including average energy) of a TPS 
very efficiently. More specifically, the relative error e in 
a TERG calculation scales with the computation time T 
like e - e-™"'''-^"^)' for 2D gapped systems. [13] (For 
comparison, the error in a variational Monte Carlo calcu- 
lation scales like e ~ e~ 2 '" -^ , if there is no sign problem) . 

In this paper, we apply the TERG approach to a few 
models with string-net condensation. We calculate the 
phase diagram of these systems and study phase tran- 
sitions from string-net condensed states to states with- 
out string-net condensation. These transitions are exam- 
ples of continuous phase transitions between phases with 
different topological orders but the same symmetry. [2-7] 
As such, they are beyond the Landau symmetry break- 
ing paradigm. Thus the TERG approach is capable of 
describing phases and phase transitions that cannot be 
described by Landau's symmetry breaking theory. 

Z2 gauge model: The first system that we study is a 
spin- 1/2 system where the spins live on links of a square 
lattice. The Hamiltonian is given by 
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Here IIigd '^i ^^ ^^^ product of the four af around a 



square p and ^ sums over all the squares. Yii^v ^i ^^ 
the product of the four ai around a vertex v and J2v 
sums over all the vertices. J^i sums over all links. We 



will assume that U — 00 and study the quantum phases 
of the above system as we change g and J. We will 
assume J > and g > 0. 

When J = 0, eqn. (1) is exactly soluble. [14] To under- 
stand the exact ground state in the string language, [15] 
we interpret the a^ = —1 and a^ = 1 states on a single 
link as the presence or absence of a string. The appro- 
priate low energy Hilbert space in large U limit is made 
of closed string states that satisfy Hie^ '^i = 1 ^t every 
vertex. The ground state is simply an equal weight super- 
position of all closed string states I^Za) = Sxciosod l"'^)' 
which is called a string-net condensed state. 

When g — 0, the ground state is the spin polarized 
state with no down spins and no closed strings. The 
above two states have the same symmetry. But due to the 
non-trivial topological order in the string-net condensed 
state, the two states belong to two different quantum 
phases. We would like to use the TERG approach to 
study the phase transition between the above two states 
with different topological orders. 

We would like to mention that the low energy effec- 
tive theory of eqn. (1) is Z2 gauge theory.[14, 16, 17] The 
transition between the string-net condensed and non con- 
densed phases is nothing but the transition between the 
deconfined and confined phases of Z2 gauge theory. 

One way to study such a phase transition is to intro- 
duce a variational wave function 
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where Lx is the number of links on the string X. When 
w — 1, ]^^,) becomes the string condensed state J^Za)- 
When w = 0, \'^w) is the state with all spins in up direc- 
tion which does not contain any strings. We see that w 
describes the string tension. ]^^=i) is the ground state 
of eqn. (1) when J = and ]^^=o) is the ground state 
of eqn. (1) when g — 0. 

Since ]\E'u,=i) and l^m=o) have the same symmetry, 
one might expect that as we change g/J, \'^w=i) would 
change into l^t^j^o) smoothly and the ground state energy 
of eqn. (1) would be a smooth function of g/J, implying 
that there was no quantum phase transition. In fact, we 
will see below that the ground state energy of eqn. (1) 
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FIG. 1: (a) The spin-1/2 model eqn. (1) on links of a square 
lattice. The dots represent the physical spin states which are 
labeled by m = 0, 1. The above graph can also be viewed as 
a tensor-network where each dot represents a rank-3 tensor g 
and each vertex represents a rank-4 tensor T. The two legs of 
a dot represent the a and (3 indices in the rank-3 tensor g^p. 
The four legs of a vertex represent the four internal indices 
in the rank-4 tensor Tap-yx- The indices on the connected 
links are summed over which define the tensor trace tTr. (b) 
A tensor-network of the double-line tensors, where each dot 
represents a double-line tensor g and each vertex represents 
a double-line tensor T. The four legs of a dot represent the 
ai,2 and /3i,2 indices in the tensor ff™j^a2;/3i/32- "^^^ eight legs 
of a vertex represent the internal indices in the rank-4 tensor 
rQiQ2;/3i/32;7i72:Ai^2- The indices on the connected links are 
summed over which define the tensor-trace tTr. 



is not a smooth function of g/J indicating that there is 
quantum phase transition at a critical value {g/ J)c- 

In order to calculate the energy expectation values in 
these states (and also to pave the way for generalizations) 
it is convenient to write the trial wave function |4',„) as 
a tensor product state: 

l*22>= Y. tTr[®„r®,g'"']|mi,m2,...), (3) 



where m; = 0, 1 labels the up-spin state and the down- 
spin state on link-L To define the tensor-trace (tTr), one 
can introduce a graphical representation of the tensors 
(see Fig. la). Then tTr means summing over all indices 
on the connected links of tensor-network. The Z2 string- 
net condensed ground state that we discussed above is 
given by the following choice of tensors: 
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1, if a + P + J + 6 even 
0, if a + I3 + J + 5 odd 



ffoo = 1, 



5i\ 



others — 0, 
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with internal indices like a running over 0, 1. The rank- 
3 tensor g behaves like a projector which essentially set 
the internal index equal to the physical index so that 
a = 1 represents a string and a = represents no string. 
The meaning of the tensor Tap-fS is also clear, it just 
enforces the closed string constraint, only allowing an 
even number of strings to meet at a vertex. 
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FIG. 2: (Color online) The black squares are average energies 
of the Z2 model eqn. (1) for the single-line tensor network 
eqn. (4) and eqn. (5). The red dots are average energies for 
the double- line tensor network eqn. (8) and eqn. (7). 



Once we have expressed the trial wave function as a 
TPS, we can use the TERG mcthod[13] to calculate the 
average energy in a very efficient way. [20] The resulting 
average energy as a function of g/J is plotted in Fig. 2. 
From the discontinuity in the slope, we see that there is 
a first order phase transition at g/J « 2.3 between the 
two states with and without string-net condensation. 

How good is this result? On a quantitative level, it is 
not very good: the phase transition is known to occur 
at {g/J)c ~ 3.044. [18] However, this not surprising since 
we used the simplest possible variational wave function. 
We expect the estimate for {g/J)c to improve when we 
increase the number of variational parameters - for ex- 
ample, by considering more general tensors g, T . 

A more serious problem is that the result is wrong 
on a qualitative level: the phase transition is known to 
be second order, not first order. This problem cannot 
be overcome by blindly generalizing the tensors g, T. 
Instead, we have to choose these tensors in a special way. 
To understand the basic issue, let us consider another set 
of variational tensors. In this scheme, the internal indices 
for the T-tensors and g-tensors still run from to 1, but 
each leg now has two internal indices: [21] 

J^afUjX ^ -'Q;ia2;/3i/32;7i72;AiA2 

= ^ai/3i7iAi'^a27l'^72/3i'5/32Ai<5A2ai (6) 

with ai,a2,;9i,;52,7i,72, Ai, A2 = 0, 1 and 

5?iai = 5^0,00 = 1> ffioao =501,01 = 1' others = 0. (7) 

In such constructions, our tensors have a double-line 
structure (see Fig. lb). Again, g™ are projectors that 
relate the internal indices with physical indices. Here, on 
each leg of g-tensor and T-tensor, the double line with 
the same value is projected to the spin-up state and the 
double line with different values is projected to the spin- 
down state. 







FIG. 3: (Color online) (Bp) versus g/J. The single-line varia- 
tional wave function shows a jump in {Bp}, which indicates a 
first order phase transition around g/J = 2.3 (black squares). 
The double-line variational wave functions has no jump in 
(Bp) but the discontinuity in the derivative indicates a sec- 
ond order phase transition around g/J — 3.2 (red dots). 



To maintain the 90 degree rotational symmetry, we 
choose r° to have a form (assuming Tim ~ !)• 



x(l) 
x{2) 
x(3) 
x(4) 
x(5) 



t; 




0000 

rpO _ rpO 

^0111 ~ ^1011 



7^0 
-'llOl 



7^0 
-'lllO 



rpO rpQ rj^O rj^O 

^1000 — ^0100 ^ -'ooio — -'oooi 
^1100 — ^0011 ^ -'oiio — -'lOOl 

rpO _ rpO 

^1010 — ^0101 



(8) 



We note that for such a choice of r° and 5™, the trial 
wave function contain only closed string states. 

Using the TERG approach to minimize the average en- 
ergy, we find the variational ground state energy which 
is plotted in Fig. 2. We find that there is a phase transi- 
tion between the two phases with and without string-net 
condensation. But now the phase transition is a second 
order phase transition at g/J « 3.2 (see Fig. 3). Note 
that this result is better than our previous result both 
quantitatively and qualitatively. 

The quantitative improvement is perhaps not surpris- 
ing since we are using more variational parameters. A 
more important issue is that the double-line mean field 
theory correctly predicts a second order phase transition, 
while the single line mean field theory did not. Why is 
this? 

Note that there is a Z2 redundancy in the double-line 
tensors (like the gauge redundancy in gauge theory). As 
we exchange values of and 1 for all the internal indices of 
the double-line tensors, we induce a Z2 transformation on 
those double-line tensors: {T,g™) -^ {T,g"^). However, 
such a Z2 transformation does not change the physical 
wave function: tTr[(8)„r®ig"'] = tTr[(g)„T(g)/g"']. Thus 
{T,g™) and {T,g"'-) are two labels that label the same 
physical state. 
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FIG. 4: The double-semion model on the honeycomb lattice. 
The ground state wave function (12) has a TPS representation 
given by the above tensor-network. Note that T and g has a 
double-line structure. The vertices form a honeycomb lattice 
which can be divided into A-sublattice and B-sublattice. 



The variational approach used here is similar to cal- 
culating an average in a local classical statistical system. 
The presence of a Z2 symmetry allows a classical sys- 
tem to have a Z2 symmetry breaking transition which 
is a second order phase transition. This is the reason 
why the double-line tensors are capable of producing a 
second order phase transition. In contrast, for the single- 
line tensors, the corresponding classical system does not 
have any symmetries, and as a result, it cannot describe 
a second order transition. 

We would like to mention that there is a duality trans- 
formation that relate the 2D Z2 gauge theory to trans- 
verse field Ising model. [19] Such a duality mapping al- 
lows us to relate the phase transition between the decon- 
fined and confined phases of the Z2 gauge theory to the 
spin ordered and disordered transition in the transverse 
field Ising model. This is how we know that the transi- 
tion between the string-net condensed and non condensed 
phases is a second order phase transition and that it oc- 
curs at critical couphng g/ J k, 3.044. In fact, the double- 
line tensors exactly realize the duality mapping between 
the 2D Z2 gauge theory and transverse field Ising model. 
From the structure of the double tensors in eqn. (7) and 
eqn. (8), we see that each square loops in Fig. lb car- 
ries the same value of internal indices, which correspond 
to the value of a dual spin (located at the center of the 
square) in the dual Ising model. The string formed by 
the down-spins on the links correspond to a domain wall 
in the dual Ising model. 

Double-semion model: Next we consider a more com- 
plicated model where spins are located on the links of a 
honeycomb lattice (see Fig. 4): 



Ices of / 



5E( n ^j)( n i^)' (9) 

p edges of p legs of p 



where i labels the links, I labels the vertices and p labels 
hexagons. Again we consider U — oo limit. When ,7 = 0, 
the above model is exactly soluble and the exact ground 
state is given by[ll] |*,em) = Ex(-)'^^^l^)' where ^^ 
sums over all the closed string configurations and 1{X) is 
number of closed loops in X. The ends of string in such a 
state have the semion statistics. When g = 0, the model 
is also exactly soluble and the spins all point up {ie no 
strings) in the ground state. 

To study the phase transition between the above two 
states, again we choose the double-line tensors to con- 
struct the trial wave function (see Fig. 4) . The T-tcnsors 
in the vertices are given schematically by 



sublattice A : 
sublattice B : 



-^ aa';/3/3';77' 






where each internal index represented by one of the dou- 
ble lines runs over 0, 1. The tensor T and T^ is given 

by 
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(11) 



The (jf-tcnsors on the links are given by eqn. (7). The trial 
wave function is obtained by summing over all the inter- 
nal indices on the connected links in the tensor network 
(see Fig. 4): 
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5^5"']|"^l,m2,...). 



(12) 



Again, the physical indices and the internal indices have 
a similar relation as in the Z2 double-line tensors. When 
x(l) = a;(5) = -i, x{2) = x(6) = i and x(3) = a;(4) = 
x(7) = x(8) — 1, the tensor reproduces the right sign 
oscillations (— )'(^) essentially by counting the number 
of left and right turns made by the string. 

We used the TERG approach to find the variational 
ground states for different g/J. Then we used the TERG 
approach to calculate {Bp) for those variational ground 
states. The result is presented in Fig. 5. We see that 
there is a second order phase transition at g/J k, 5.0, 
which should correspond to the transition between the 
string-net condensed and non condensed states. This 
agrees with the Monte Carlo result where a second or- 
der phase transition appears at {g/J)c ~ 4.768. [18] (Note 
that in the infinite U limit, the above model is equivalent 
to the usual Z2 gauge model on honeycomb lattice, which 
is dual to the transverse Ising model on triangle lattice.) 

Conclusion: We have seen that the TERG approach is 
an effective way to study topological phases and topolog- 
ical phase transitions, but one needs to choose the varia- 
tional tensors carefully. An important question is how to 




FIG. 5: (-Bp) versus g/J, the discontinuity of the derivative 
around g/J = 5.0 indicates a second order confinement-de- 
confinement phase transition. 



choose the tensors in more general situations. One hint 
is that all the string-net states constructed in Ref. 11 
can be expressed naturally in terms of a generalization 
of the double-line tensor network, which involves triple 
line tensors. [12] This triple-line tensor-network may cor- 
respond to the dual representation of the string-net states 
and may correspond to a suitable choice for the varia- 
tional TERG approach. This may lead to a systematic 
variational approach for topological phases and topolog- 
ical phase transitions. 
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